Neutron scattering study of the field-induced soliton lattice in CuGeC>3 
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CuGeOa undergoes a transition from a spin-Peierls phase to an incommensurate phase at a critical 
field of H c ta 12.5 T. In the high-field phase a lattice of solitons forms, with both structural and 
magnetic components, and these have been studied using neutron scattering techniques. Our results 
provide direct evidence for a long-ranged magnetic soliton structure which has both staggered and 
uniform magnetizations, and with amplitudes that are broadly in accord with theoretical estimates. 
The magnetic soliton width, F, and the field dependence of the incommensurability, Sksp, are found 
to agree well with theoretical predictions. 



A spin ~ antiferromagnetic (AF) chain in a deformable 
3D lattice may, at a characteristic temperature T SP , un- 
dergo a spin-Pcicrls transition into a correlated dimer- 
ized groundstate characterized by a total spin Stat = 0- 
The energy cost of deforming the lattice is more than 
compensated for by the reduction in magnetic energy, 
and a gap opens in the magnetic excitation spectrum 
which reduces the influence of zero-point fluctuations [0 . 
This intriguing magneto-elastic phenomenon gives rise 
to several interesting properties. The groundstate is pre- 
dicted to be a non-magnetic singlet separated from an ex- 
cited triplet (5tot = 1) by a gap of energy A = 1.76 T SP . 
The opening of this gap produces anomalies in the tem- 
perature dependence of physical quantities like the mag- 
netic susceptibility, specific heat, etc. Dimerization also 
results in a doubling of the unit cell along the chain di- 
rection, which is reflected in diffraction experiments by 
the appearance of satellite reflections at commensurate, 
half-order positions. 

Application of a magnetic field reduces the stability of 
the spin-Peierls (SP) state against a state with finite mag- 
netization and, at a critical field H c ~ 0.84Ao/g/is, the 
system enters an incommensurate high- field phase |],[j). 
Above H c , the development of a uniform magnetization 
can be visualized as the breaking of a dimer bond, thus 
allowing two spins to become parallel to the magnetic 
field. The two S z = \ states can be regarded as a pair of 
domain walls with respect to the dimer order. The trans- 
verse part of the AF interaction delocalizes the domain 
walls into solitons of finite width T, and the longitudi- 
nal part causes the two domain walls to repel each other, 
thus forming an equally spaced lattice of solitons. The 
resulting soliton has both ferromagnetic and AF compo- 
nents and, in addition, the lattice distortion also becomes 
incommensurate. The average magnetization is propor- 
tional to the number of solitons, and inversely propor- 
tional to the separation L. 

Theoretically, the coupling of the spin chains to the 
lattice is usually treated within the adiabatic approxima- 



tion. A continuum model has been developed by apply- 
ing bosonization techniques to the Jordan- Wigner trans- 
formed Hamiltonian ||). This field-theoretical approach 
provides an analytic description of the soliton structure 
in terms of the soliton spacing, L/c = l/5k SP , soliton 
width, r, and the amplitudes of the uniform, m u , and 
staggered, m s , magnetizations [||-|]]. The formalism has 
been extended to include next-nearest neighbour (nnn) 
and inter-chain interactions J^,|^] , but it is somewhat un- 
clear to what extent the results are affected by the various 
approximations involved. Monte Carlo |[o| and density 
matrix renormalization group (DMRG) calculations have 
also been performed 0,0] . While the soliton shape 
matches the field theoretical solution, the DMRG cal- 
culations show that the structural and magnetic soliton 
widths are in general different to each other and field 
dependent. 

As the high-field soliton lattice has both a structural 
and magnetic component it is desirable to use a probe 
that is sensitive to both. In this letter we report the re- 
sults of an investigation of the structure of the high-field 
phase of CuGc03 using neutron scattering techniques. 

Experimental investigations of the SP transition were 
initiated by the work of Bray et at on the organic com- 
pound TTFCuS 4 C 4 (CF 3 )4 || . Interest in SP systems 
intensified with the discovery by Hase et al. |1| that 
CuGeOa undergoes an SP transition at T SP = 14 K. 
CuGeOa crystallises in the orthorhombic Pbmm struc- 
ture with room-temperature lattice parameters of a — 
4.80 A, b = 8.47 A and c = 2.94 A, and the S = \ Cu 2+ 
moments form chains along the c-axis. The field theory 
predictions depend on the spin gap Ao = 2 meV and 
spin-wave velocity v s = 16.6 meV, which have been de- 
termined using inelastic neutron scattering |l4|-|l6|]. The 
inter-chain couplings J a = —0.08 meV and Jb = 0.62 
meV indicate that the system is in fact only quasi- ID 
fL7j| . Along the chain, both an alternating exchange cou- 
pling J± = (1 ± 0.014) x 14 meV and a significant nnn 
coupling J' c ~ 0.35J C are needed to model the suscepti- 
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bility, saturation field and excitation spectrum 

In zero-field the expected doubling of the unit cell 
along the c direction below T SP has been confirmed by 
electron ]is|] , X-ray |^(| and neutron pl[ , p2t diffraction 
experiments. The strongest satellite reflections appear at 
(5, k, 5) with h, k and I odd, which indicates that along a 
and b adjacent chains dimerize in anti-phase. Hirota et al. 
pl[ and Braden et al. Q have solved the detailed struc- 
ture of the SP phase, which is characterized by the prop- 
agation wavevector fc SP = (§,0, 5) and three displace- 
ment parameters: u° u / c = 0.00192, u° 2 /a = 0.00198, 
M Q2/& = 0.00077. 

Several experiments have been reported on the high- 
field phase, including X-ray scattering f^jj , thermal ex- 
pansion §J, FIR spectroscopy ||, NMR ||, ESR and 
magnetization f27j , p8|| . The predicted incommensurabil- 
ity of the structure was observed in an X-ray experi- 
ment [^3| as a splitting of the satellite reflection into 
(5, 1, | ± 6k SF ), and the soliton width was estimated to 
be r = (13.6 ± 0.3)c. Using copper NMR, Horvatic et 
al. |2(| have determined the staggered and uniform mag- 
netizations to be m s = 0.026 and m u = 0.023, respec- 
tively at 14.5 T. (By definition a spin of 5 corresponds 
to vci — i.) Magnetization measurements give a value of 
0.018 for m u . The theoretical values are m s = 0.11 <-» 
0.14 and m u =0.023 Thus although there is good 

agreement between theory and experiment for to u , the 
staggered moment observed in NMR is significantly lower 
than the theoretical value, a point we will return to later. 

Although neutron scattering is an obvious way to at- 
tempt to study the magnetic soliton lattice, it has until 
recently not been possible to apply this technique above 
12 T. In CuSi.03Ge.97O3, where H c is reduced to 11.7 
T, Grenier et al. |2S| ] have measured the intensity of five 
satellite reflections upon entering the high-field phase. 
No clear evidence of any magnetic scattering was found. 
The soliton width at 12 T was deduced to be 13.5c, in 
agreement with the X-ray data on pure CuGe03. 

Here, we report on a neutron scattering study of the 
high-field phase in pure CuGeOs, performed at HMI, 
Berlin, using a 14.5 T vertical field cryomagnet with a 
base temperature of 1.6 K. The single crystal of dimen- 
sions 3 x 7 x 16 mm 3 was grown by the floating zone 
method. In total three different instruments were used, 
each chosen to reveal a particular aspect of the phase 
transition. Due to the slightly anisotropic g-tensor [ j30| , 
the critical field depends on the orientation of the applied 
magnetic field. 

The field dependence of the incommensurability was 
investigated using the thermal triple-axis spectrometer 
El with a collimation of 20'-40'-20'-80', and a wavelength 
of A = 2.425 A. By performing scans along c* through five 
satellite reflections (5,1,5), (5,3, 5), (5,3,5), (5,1,5) 
and (5,3,5), we observed the expected splitting. 
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FIG. 1. The measured splitting Sk SP of 1, |) (□) and 
(5,3, |) (o), compared to X-ray data ^] (ic), and to the 
theories of Cross and Fisher |2j (dashed line), and Buzdin 
et al. Q (solid line with Ao/v a = 0.13). For clarity, only a 
representative selection of the X-ray data is shown, and all of 
the data has been scaled to the g factor corresponding to the 
(h, 2h, I) orientation. The insert shows the soliton width V as 
obtained from the relative intensity of the 3rd harmonics. 

As shown in Fig. |[ the field dependence of 5k SP is con- 
sistent with the X-ray results in the region of overlap, 
but our data extend to much higher fields. Field theory 
predicts that the incommensurability of the soliton 
lattice is given by 1 /Sk SP — ^ In ^ Hl H . Using the up- 
per limit value of 0.13 for the experimentally determined 
ratio Ao/vs = 0.12 ± 0.01, we obtain perfect agreement 
with our data. At the highest fields, the linear field de- 
pendence Sk SP = g/iBH/ (2irv s ) is approached [0. 

One striking feature of this data set is that there was a 
tremendous increase in intensity of the low Q satellites. 
This appears difficult to explain on the basis of a purely 
structural distortion, for which the intensity increases in 
proportion to Q 2 . One possible explanation could be 
that the magnetic soliton lattice has the same propaga- 
tion vector as the lattice soliton. At first sight this seems 
unlikely, as the minimum energy gap in the spin-wave 
spectrum is located at (0, 1, 5), which is also the propa- 
gation vector of the AF structure found in doped systems 
pjj . It should be noted, however, that the magnetic cou- 
pling in the a direction is weak, of order 1 K, and could be 
overcome by stronger lattice forces. To exclude the pos- 
sibility that there is magnetic scattering from the stag- 
gered component at wavevectors other than (5, k, |), an 
exhaustive search for additional satellite reflections that 
may be associated with long range magnetic order was 
conducted. From this an upper limit can be put on any 
Fourier component at (0,1,5), (0,3,4) and (0,0, 5) of 
0.005/xs at 14.5 T. 

To investigate whether in fact there is a magnetic con- 
tribution to the low-Q satellites, a second experiment 
was undertaken. The two-axis spectrometer E4 was used 
with a collimation of 40'-40'-40' and A = 1.2205 A, which 
made 14 satellite reflections accessible in the (h, 2h, £)- 
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plane (see Table ||). Rocking curves of each reflection 
were collected at 2 K in zero field and 14.5 T. An addi- 
tional data set was taken at 20 K (above T SP = 14.1 K) 
and T to determine the A/2 contamination. The widths 
of all measured peaks were in complete agreement with 
an analytic resolution calculation, and this was used to 
correct the integrated intensities of the measured satellite 
peaks. The overall scale factor that converts the satellite 
intensities to structure factors can be found by accurately 
determining the structure factors of the main Bragg 
peaks. However, this process is often problematic, as 
it involves additional corrections for extinction, absorp- 
tion, and in our case A/2 contamination. To avoid these 
problems, we chose instead an overall scale factor such 
that the sum of the displacement parameters obtained 
when fitting our zero-field data matched the sum of the 
parameters obtained in the detailed study by Braden 
et al. |p2fl . The relative size of our zero-field displace- 
ment parameters uj u /c = 0.00204 ± 0.00014, u° 2 /a = 
0.00182 ± 0.00010 and u° 2 /6 = 0.00079 ± 0.00009 agree 
well with the values reported by Braden et al. 
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FIG. 2. Rocking curve through (|,3, \ + Sk SP ) at 14.5 T, 

2 K. The finite instrument resolution picks up both satellites 

(|,3, !±<5fc SP ), and the A/2 background at (|,3, |). The solid 

line shows a fit of three resolution convoluted peaks. The 

insert displays the reciprocal space trajectory of the scan and 

the resolution ellipse. The position of each peak is marked 

with a cross. 



The same scale factor was used for the data taken in 
the high-field phase at 14.5 T (Table @). The resolution 
obtainable with the short wavelength was insufficient to 
completely separate the two satellites (§, ft, |±5fc SP ) and 
the A/2 background at all (|, ft, |). This is illustrated in 
the scan through (|, 3, ^+6k SP ) shown in Fig. [| The data 
were fitted to three peaks convoluted with the calculated 
resolution. The A/2 contribution was fixed to the value 
determined at 20 K and T, and the two satellites were 
given equal amplitude. With the satellite amplitude as 
the only fitting parameter, all fourteen reflections gave 
excellent fits as shown by the solid line in Fig. [|. The 
resulting structure factors are listed in Table |. 
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TABLE I. Measured and fitted structure factors in units of 
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As anticipated, it was not possible to obtain satisfac- 
tory fits of the high-field structure factors by assum- 
ing a structural distortion alone. Instead, a magnetic 
Fourier component of variable size and direction with 
the magnetic form-factor of Cu 2+ was included in the 
fit. The resulting displacement parameters u^ u /c — 
0.00158 ± 0.00026, u° 2 /a = 0.00196 ± 0.00036 and 
u° 2 /b = 0.00061 ± 0.00025 are not significantly differ- 
ent from the zero field parameters. The direction of 
the magnetic moments was found to be parallel/anti- 
parallel to the field, in agreement with theoretical ex- 
pectations. The Fourier component was determined to 
be fi(q = Sk SP ) = 0.098 ± 0.003/i B at 14.5 T. 
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FIG. 3. Scans along [00/] through (|, 1, |) at 1.6 K and 
T (squares), and 13 T (filled circles). The resolution limited 
peaks were fitted to Gaussian line shapes. The higher har- 
monics n5k SP are marked by arrows. The insert shows scans 
through (0,0,1 + 2Sk SP ) at 13 T (filled circles) and 14.5 T 
(diamonds). 



Having established that the (5, 1, \ ± 8k SF ) reflection 
is almost completely magnetic in origin, we performed 
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a high-resolution experiment in order to investigate the 
detailed structure of the magnetic soliton lattice. Using 
the cold triple axis spectrometer V2 with A = 5.236 A 
and 20' collimation gave a resolution of 0.009 A -1 . Fig. 
U shows scans through (|, 1, |) along [00/] at T and at 
13 T. Along with the dramatic increase in intensity, the 
13 T scan reveals both third and fifth harmonics, which 
reflect the non-sinusoidal shape of the magnetic solitons. 

From field theory, an analytical solution for the magne- 
tization {S z (l)) = m s (-l)'cn(7c/rfc,fc) + m u dn(Zc/rfc,fc) 

at site I with amplitudes m s = « / w^J- and m — A 
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is given in terms of the Jacobian elliptic functions cn 
and dn of modulus k, where L/T = 4kK(k) and K(k) 
is the complete elliptic integral of first kind ||-||]. The 
nth Fourier components of the staggered and uniform 

n/2 

parts of the magnetic structure are given by j^jrs i +K n 

n/2 

(n odd) and K ^ i +K n {n > even) respectively, where 

K = e -^K{y/l=W)/K{k)_ p rom the ratio K 2jQ _ K + K 2j2 

between the intensities of the third and the first har- 
monic, we determine the modulus k and the soliton width 
T/c = l/(45k SP kK(k)), which is shown in the insert of 
Fig. [l]. Above the transition, T/c quickly decreases from 
11.5 to a value of 9, followed by a slight increase with 
field, a feature which is also apparent in DMRG calcu- 
lations p^ |. Except for the rapid change just above H c , 
our data are in fair agreement with the constant value 
T/c = v s /A = 8.3 predicted by field theory, but dif- 
fers from the values T/c ~ 13.5 obtained from X-ray 
and neutron scattering on the Si doped system. We be- 
lieve that this is because the latter experiments measure 
the width of the structural soliton, which in the pres- 
ence of nnn coupling is larger than the magnetic soliton 
width Jl2| . Re-analysing the X-ray data, by allowing for a 
field dependent soliton width, leads to a similar decrease 
T/c = 15.1 — > 12.9 with increasing field. 

Knowing the precise soliton shape, the amplitude 
of the staggered moment can be determined from the 
Fourier component: m s = kK{k) l±g ^"= 5k ^ = 0.097 ± 
0.003, where g = 2.19 for B _L (h, 2h, I) @. This value is 
in reasonable agreement with theoretical estimates cited 
earlier, but is considerably larger than the value derived 
from NMR. If, as has been suggested, the NMR value 
is reduced by quantum fluctuations in the groundstate 
p6| , ^2[ , then these would also reduce the value of m s ob- 
served in our neutron scattering experiments. The NMR 
and neutron scattering data can be reconciled if there are 
fluctuations fast compared to the time scale probed by 
NMR (~ 10" 4 THz) but slow compared to that probed 
by the neutron scattering (~THz). 

The uniform part of the magnetic structure should give 
rise to even harmonics around integer positions in recip- 
rocal space. Indeed, second-order harmonics were found 
near the (0,0,1) Bragg reflection, as shown in the insert of 
Fig. ||. To extract the absolute structure factor, we scale 



the integrated intensity to the (|, 1, \ ± 5k SP ) reflection 
measured in the same configuration, correcting for the 
different resolution volumes. The resulting Fourier com- 
ponent of the magnetic structure is fi(0, 0, 1 + 2Sk SF ) = 
0.0046^b at 14.5 T. The amplitude of the uniform com- 

, • • i K(k) 1 + k, 2 «(g=2(5fc SP ) 

ponent is then given by m u = ^ * g — = 
0.019 ± 0.003, about 25% lower than the field theory and 
NMR results. It should be noted that although the indi- 
vidual values of m s and m u are sensitive to the details of 
any normalisation procedure, the ratio m s /m u — 5.1±0.8 
is not, and is seen to be in good agreement with the the- 
oretical ratio m s /m u = 5.4.. 

In summary, neutron scattering has been used to de- 
termine the nature of the soliton lattice in the high-field 
phase of the spin-Peierls compound CuGeC>3. The re- 
sults provide a detailed description of the soliton lattice 
in terms of incommensurability, Sk SP , soliton width, F, 
and the amplitudes m s and m u of both the staggered 
and uniform parts of the magnetization. 
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